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Abstract 

We consider the dynamics of small perturbations of stable two-frequency 
quasiperiodic orbits on an attracting torus in the quasiperiodically forced 
Henon map. Such dynamics consists in an exponential decay of the radial 
component and in a complex behaviour of the angle component. This be- 
haviour may be two- or three-frequency quasiperiodicity, or it may be irreg- 
ular. In the latter case a graphic image of the dynamics of the perturbation 
angle is a fractal object, namely a strange nonchaotic attractor, which ap- 
pears in auxiliary map for the angle component. Therefore, we claim that 
stable trajectories may approach the attracting torus either in a regular or in 
an irregular way. We show that the transition from quasiperiodic dynamics 
to chaos in the model system is preceded by the appearance of an irregu- 
lar behaviour in the approach of the perturbed quasiperiodic trajectories to 
the smooth attracting torus. We also demonstrate a link between the evolu- 
tion operator of the perturbation angle and a quasiperiodically forced circle 
mapping of a special form and with a Harper equation with quasiperiodic 
potential. 

PACS numbers: 05. 45. -a 
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I. INTRODUCTION 



Prom the beginning of 1980s, much attention has been paid to studies of quasiperiodically 
forced systems. Indeed, systems with controllable ratios of incommensurate frequencies 
represent convenient models for the analysis of bifurcations of quasiperiodic regimes and 
of the mechanisms of transition from quasiperiodicity to chaos. Attention to this class of 
systems emerged also due to strange nonchaotic attractors-*^ (SNAs), which can generically 
occur there. A large number of works was devoted to the investigation of mechanisms 
of regular bifurcations [12-19], irregular dynamical transitions [15-25] and attractor crises 
[26-33] in quasiperiodically forced systems of different nature. 

In most cases, interest in the dynamical systems is focused on the analysis of station- 
ary dynamical regimes and their transformations. Researchers deliberately leave aside the 
dynamics of transient processes. In our humble opinion, such an approach is not wholly jus- 
tified, since transient trajectories may visit wide regions of the phase space before drawing 
near the attractor. As a consequence, transient processes may contain extra information 
about the structure of invariant sets in the phase space of any dynamical system. Moreover, 
a change of the character of the transient process under variation of the controlling param- 
eters of the system may precede transformations of the stationary dynamical regimes. In 
particular, a change of the character of the transient process may signal possible bifurcations 
and crises of attractors under further variation of the parameters of the system^. 



^Strange nonchaotic attractors were first described in Ref. [1]. The term "strange" characterizes a 
fractal-like geometrical structure of the attractor, while the term "nonchaotic" suggests absence of 
exponential instability of the trajectories on the attractor. Appearing on the border of quasiperi- 
odicity and chaos, attractors of this type possess mixed properties of the both types of dynamics. 
For more details on properties of the SNA see Refs. [2-11]. 

^Refer to the node and focus stable fixed points of a 2D map. Linear analysis of evolution 
of small perturbations near such fixed points reveals different character of approach of transient 
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Here, we proceed with the investigation of the bifurcations and dynamical transitions in 
quasiperiodically forced systems, although we focus on dynamics of transient processes and 
their transformations. Our main goal is to show that a regular attractor may have an irregu- 
lar transient process in its vicinity, and a complication of the transient may signal imminent 
destruction of the regular attractor under variation of the parameters of the system. 

In the present work we consider the dynamics of small perturbations of stable quasiperi- 
odic orbits with two incommensurate frequencies on an attracting torus. As a model system, 
we consider the quasiperiodically forced Henon map, which can be interpreted as a Poincare 
map of a hypothetical nonlinear oscillator driven by external biharmonic signal with irra- 
tional frequency ratio. A smooth closed invariant curve in the 3D phase space of the model 
map corresponds to a Poincare section of a 2D torus in the 4D phase space of such an 
oscillator. Therefore, we shall refer to such a smooth closed invariant curve as a "torus". 
Evolution of a small perturbation of a trajectory on a stable torus consists of exponen- 
tially decaying radial component of the perturbation (what corresponds to the approach of 
the perturbed trajectory to the attracting invariant curve) and complex behaviour of the 
angle component of the perturbation, which characterizes a rotation of the perturbation 
vector around the invariant curve. We show that the dynamics of the angle component of 
the perturbation may have the character of a two- or three-frequency quasiperiodicity, or 
it may be irregular. In the first two cases a graphic image of such dynamics represents a 
two- or three-frequency torus, while in the latter case this image corresponds to a strange 
nonchaotic attractor. Note, that the attractor of the model system remains a smooth torus, 
while irregular transient process appears in its vicinity. 

We analyse the global structure of the parameter space of the model system in the 



trajectories to different fixed points. For the stable fixed points of these two types, a small variation 
of parameters may give rise to different local bifurcations: doubling, pitchfork, inverse saddle-node 
etc. for nodal fixed points, and Neimark bifurcation only for focus fixed points. 
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region of complex dynamical transitions consisting in the destruction of the smooth torus 
accompanied by the birth of a strange nonchaotic attr actor or a divergence of trajectories. 
We demonstrate that the appearance of an irregular (in terms of angle) character of approach 
of stable trajectories to the attracting torus plays the role of a "precursor" of such transitions. 
Namely, the birth of a SNA via torus gradual fractahzation [20,21], intermittency [23,24] 
and collision with a saddle torus [22], as well as the appearance of trajectory divergence 
via the collision of an attracting torus with a fractal basin boundary [26,27], is preceded by 
the appearance of irregular dynamics of the perturbation angle. Note that the attractor of 
the system remains smooth while perturbed trajectories in its vicinity already demonstrate 
"strange" behaviour. On the other hand, regular torus bifurcations such as torus-doubling 
are accompanied by regular behaviour of the perturbation angle both before and after the 
bifurcation. We also discuss aspects of a link between the perturbation angle dynamics and 
a special circle mapping and with the Harper equation with a quasiperiodic potential. 

The paper is organized as follows. In Sec. II we briefly observe the types of attractors 
and dynamical transitions in the parameter space of the model system. In Sec. Ill we obtain 
the equation of the evolution of the perturbation angle in explicit form, and discuss its 
connection with the Harper equation (Sec. IV). In Sec. V we analyse numerically different 
kinds of the perturbation angle dynamics and demonstrate their relationship with the types 
of dynamical transitions in the parameter space of the model system. 

II. BIFURCATIONS AND DYNAMICAL TRANSITIONS IN THE 
QUASIPERIODICALLY FORCED HENON MAP 

Our model system has the form: 

Xn+l = f{Xn, Vn: On) = a - xl - by^ + S COs6l„, 

Un+l = g{Xn, Vn, On) = Xn, (l) 

On+1 = On + 2t^uj (mod 27r), 
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where the quasiperiodic force frequency parameter is chosen equal to the inverse golden 
mean value: to — — l)/2. 

The examples of typical attractors of the map (1) are presented in Fig. 1. At a = 1.09, 
e — 0.2 (here and hereafter we set b — 0.3) the attractor is a smooth closed invariant curve 
("torus" T, Fig. 1(a)). With an increase of the nonlinearity parameter a, the torus T may 
bifurcate into a "double torus" 2T, which consists of a pair of smooth closed curves mapping 
into each other under iteration of the map, as shown in Fig. 1(b) where a — 1.445, e — 0.095. 
(If £ is small enough, one can also observe a second torus-doubling bifurcation, which gives 
four smooth closed curves. Note, that the number of torus-doubling bifurcations on the route 
to chaos depends upon the quasiperiodic force amplitude e; this number may increase and 
tend to infinity as £ goes to zero [34] . The mechanism for the termination of torus-doubling 
cascades in invertible systems is discussed in Ref. [35].) On the other hand, a different 
kind of torus-doubling bifurcation (torus "length-doubling" bifurcation [29]) can occur for 
sufficiently large values of e. This bifurcation produces a special type of torus, which is 
characterised by the existence of two wraps of a single invariant curve; we will refer it to as 
a "double-wrapped torus" T2 (see Fig. 1(c) where a — 0.728, e — 0.36). Strange nonchaotic 
and chaotic attractors of the map (1) are presented in Fig. 1(d) (a = 0.94, e — 0.312) and 
Fig. 1(e) (a = 1.6, e = 0.135), respectively. 

The phase diagram of the map (1) in the a — e parameter plane {b — 0.3) is shown 
in Fig. 2(a). The regions of different types of dynamical behaviour are shown in different 
tones. Torus T (white), double torus 2T (light-grey), double- wrapped torus T2 (light-grey). 
Regions of chaotic dynamics are shown in black. On the border between quasiperiodicity 
and chaos regions of SNA (dark-grey) exist. In the patterned area the map (1) has no 
attractor, and the trajectories go to infinity. 

In Fig. 2(b) an enlarged fragment of the phase diagram in the region of complex dy- 
namical transitions is presented. Different lines passing along the border of the region of 
quasiperiodic dynamics correspond to different dynamical transitions occurring at the exit 
from this region. The dashed line F corresponds to the birth of a SNA via gradual frac- 
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talization of the smooth torus T [20]. The dotted hne I corresponds to the intermittent 
mechanism of the birth of a SNA [23] via coUision of torus T with a saddle chaotic invariant 
set [24] . The combined hne M corresponds to the birth of a SNA via band- merging colhsion 
of a double torus 2T with a saddle parent torus [22]. A transition B from the region of 
quasiperiodicity to divergence of trajectories occurs when the torus T destroys and trans- 
forms into a chaotic transient via collision with a saddle chaotic set on the fractal basin 
boundary [26,27]. Arrows give the direction of movement in the parameter space for each 
transition to be observed. 

III. PERTURBATION ANGLE DYNAMICS AND A SPECIAL CIRCLE 

MAPPING 

Let us take a smooth attracting torus of the map (1), given by invariant curve of the 
form 

T{„,,} : {{x,y,e) e R2 X = X{a,b}{e),y = Y^a,i>}{e),e e [0,27r)}. (2) 

The pair of continuous and smooth functions {Xs^a,b}{d)i Y{a,h}{d)} can be obtained by solving 
the system of nonlinear functional equations 

X(9 + 2nu;) - a - X^(9) - bY(9) +ecos9, 

(3) 

Y{9 + 27ruj)^X{9). 

In what follows we can omit mention of the dependence of these functions upon the param- 
eters {a, 6} and write simply {X{9),Y{9)}, although we should keep in mind the link of 
{X{9),Y{9)} with the system (3). 

Take a reference quasiperiodic trajectory {{xn,yn, 9n)}n=o,i,...,oo ^ T, and consider a 
perturbed trajectory close to it starting with the same initial phase: -|- -|- 

^yn,(^n)}n=o,i,...,oo- Evolutiou of a Small perturbation 5r„ — {5xn,Syn){n — 0,1,..., oo) 
of the trajectory on torus in linear approximation is given by the following map 
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6Xn+l 



J(^n) Urn ^n) 



5Xn 

5yn 



(4) 



9n+i = + 27ra; (mod 27r) , 



where 



J(x,y,^) 



df{x,y,0) df{x,y,9) 
dx dy 

dg{x,y,e) dg(x,y,e) 
dx dy 



(5) 



Coordinates x„ and y„ in the map (4) are taken at the corresponding points of the trajectory 
on the torus T: a;„ — X{9n), Un = Y{6n). For convenience, let us introduce terms 



Tp (a \ _ df{xn,yn,6n) 
^xK^n) — 

Gx{dn) 



df(Xn,yn,d„) 



Xn=X(en},yn=Y{e„) ' ^f" 



dg(Xn,yn,dn) 
dXn 



Xn=X{en), yn=Y{en) 



dVn 



x„=X{9„), y„=Y(0„) ' 

x„=x{e„), y„=Y{en) ' 



Now the map (4) can be rewritten as 



6yn+i = Gx{9n}SXn + Gy{9n)Syn, 
9n+i — On + Sttu; (mod 27r). 



(6) 



In polar coordinates 



5x — Sr cos (p, 5y — 5r sin (f 



the map (6) becomes 



5rn+l cos (Pn+l = Fxidn)Srn COS (fin + Fy{9n)Srn Sin (fin, 
5rn+lSm(fin+l = Gx{9n)SrnC0S(fin + Gy{9n)SrnSm(fn, 

9n+i = 9n + ^TTLU (mod 2tt). 
Substituting into (7) the values 

F,{9r,) = -2X{9r,), Fy{9n) = -6, 



(7) 



Gx{9n) — 1, 



Gy{9n) — 0, 



we have 



(5rn+i COS V9„+l = -5rn[2X{0n) COS ipn + bsilKfin], 

(8) 

9n+i — 9n + Stto; (mod 27r), 
and thus, eliminating the radius variable dr, we obtain 

cot(^„+i = - [2X{9n) + 6tan93„], 

(9) 

9n+i — 9n + '^TTUJ (mod 27r). 

In order to get an explicit map for the angle variable </?, one should take into account the 
conditions 

ifn e [0, 7r/2) U [37r/2, 27r) ^ ip^+i G [0, tt), ^^^^ 
e [7r/2,37r/2) ^ e [7r,27r), 

which immediately follow from the second equation of the system (8). Then, from (9) 
and (10) one can finally obtain the explicit mapping for the evolution of the perturbation 
angle variable </?: 

ipn+l^ f{(Pn,X{9n)), ^^^^ 

9ji+i — 9n + ^TTLO (mod 27r), 



fiv,X{9)) 



where 

arctan[2X(0) + 6tan(^] + (7r/2), e [0,7r/2) U [37r/2,27r), 

arctan[2X(^) + 6tan(^] + (37r/2), ip e [7r/2, 37r/2). 

Recall that X(9) is a function of period 27r. The map (11) represents a kind of a circle 
mapping under external qiiasiperiodic forcing. Note, that attractors of the map (11) must 
be symmetric with respect to a shift tt along the (/7-axis. Therefore, the map (11) could be 
simplified to the form 

(fin+i = arctan [2X{9n) + 6tan(/7„] + (7r/2), 
9n+i — 9n + ^TTLO (mod 27r), 

by the choice of e [0, tt). 

8 



Returning to the radial coordinate 5r, one can see, that its dynamics is trivial. Indeed, 
for any small perturbation in a vicinity of the attracting torus T we immediately have: 

5rn = \J {SxnY + {5ynY ~ exp (nAi) \J{5xqY + iSyoY = 5ro exp (nAi), 

where Ai < is the largest nontrivial Lyapunov exponent, which characterizes stability of 
the torus T. 

The analysis of the dynamics of the perturbation angle (/? on complex tori of the form 
fcT (A; = 2,4, 8, . . .) in the map (1) requires trivial modification of the map (11). The last 
becomes cyclic: 

9n+i — On-\- 2'KUJ (mod 27r) , 

where n — kp + m, p is integer, m = 1,2, and {X'^'^\6)} is a set of functions 

determining a set of k smooth closed curves which the complex torus k T consists of. On the 
other hand, for analysis of the dynamics of the perturbation angle (p on the double- wrapped 
torus T2 it is convenient to redefine the phase variable d on the interval [0, 47r) and to rewrite 
the map (11) in the following form: 

9n+i — 9n-\- Stto; (mod 47r) , 
where X{2){0) is a function of the period 47r defining the double- wrapped torus T2. 

IV. A LINK TO THE HARPER EQUATION 

Let us make a change of the variable in the map (9): Un — tan (/?„. Then, one can write: 

Un+l^-[hUn + 2X{en)\-\ 

(12) 

dnj^i = 9n + 2t:uj (mod 27r). 
Then, introducing Un — b^^^Un, one obtains: 



9 



Un+l = -[Un + 2b-'/^X{en)]-\ ^ ^ 

(13) 

9n+i — On + 27ra; (mod 27r) . 
Since the function X{9) is 27r-periodic, it can be decomposed into Fourier series: 

oo 

X{0) = ao/2 + ^ A cos k{0 + (t>k). (14) 
fe=i 

Fourier coefficients {>ljt}fe=i...oo and {0fe}fe=i...oo can be obtained from numerical solution of 
the functional equation (3) . For the case of small values of the quasiperiodic force amplitude 
£ (see the map (1)), one may take into account the first terms of the Fourier series only and 
obtain 

= -K + h-^'^ao + 26-1/2^1 cos [On + 0i)]-\ 

(15) 

— 9n-\- ^TTCU (mod 27r). 
Introducing parameters E = — 6~^/^ao, A = b~^^'^Ai, the map becomes 

rt„+i = -[?/„-£; + 2A cos (6'„ + 0i)]""\ 

(16) 

6'n,+i = 6*^ + 27ra; (mod 2tt), 
which reduces to the Harper equation with a quasiperiodic potential after a standard change 
of the variables m„ = ipn-i/ipn- 

ipn+i + + 2A COS (27ra;n + 0)^„ = Etpn, (17) 

where (f) — 6q -\- is & phase shift. 

The states of the wave function in the Harper equation are known to be related to the 
dynamical regimes of the map (16). The extended states in the subcritical region (A < 
1) correspond to tree- frequency tori, while the localized states in the supercritical region 
(A > 1) are associated with SNAs in the map (16). A two-frequency torus of the map (16) 
corresponds to the non-normalizable wave function and to the gap in the energy spectrum. 

In Fig. 3 the standard phase diagram of the Harper map (16) is presented. The extended 
phase (Ext.) is shown in white, and the locahzed phase (Loc.) is denoted by gray. The 
Gap regions are shown in the light-gray tone. In order to obtain this diagram, we used 
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the following technique. Each dynamical regime of the map (16) was characterized by the 
non-trivial Lyapunov exponent u and the phase sensitivity exponent b [11]. These exponents 
take values cr = and 5 = for a three-frequency torus, cr < and 5 = for a two-frequency 
torus, and cr < with 5 > for the case of a SNA. 

The hnk from the Harper equation to the map (11) suggests that the configuration of 
the phases on the diagram (Fig. 3) can be related to the structure of sets in the parameter 
space of the quasiperiodically forced Henon map. Indeed, in the next section we will show 
that regions of different perturbation angle dynamics are organized as "tongues" of two- and 
three-frequency quasiperiodicity or of a complex behaviour corresponding to a SNA. 

V. ANGULAR DYNAMICS OF PERTURBATIONS IN THE CONTEXT OF 

TORUS BIFURCATIONS 

Now we will focus on the angle dynamics of small perturbations in a vicinity of the torus 
T. Suppose we have a trajectory starting from the initial phase on the torus and examine 
its small perturbation with an arbitrarily chosen initial angle (^o- Further evolution of the 
perturbation angle is determined by the map (11). In what follows, we choose a few sets 
of the parameters (a,£) such that a smooth torus T exists in the map (1) and consider a 
dynamics of the map (11) at the same parameter values. 

At a = 0.8, £ = 0.1 the attractor of the map (11) represents a smooth torus (p = $(6*) 
shown in Fig. 4(a). Note that tori of the map (11) may have different topology, demonstrate 
one or many wraps on the circle (p e [0, 27r) and consist of one or two segments, as shown 
in the Fig. 4(b) where a = 1.3, s = 0.5. The 2D map (11) is characterized by two Lyapunov 
exponents. A trivial Lyapunov exponent associated with the quasiperiodic variable 9 is equal 
to zero, while a nontrivial Lyapunov exponent in our cases takes the value a — —0.1170 for 
the first set of parameter values and a = —0.9755 for the second set. The angle variable ipn of 
an arbitrarily chosen perturbation on the torus T asymptotically behaves as (</?„ — $(^n)) ~ 
exp{an), tending to $(^n) as n — > oo. Thus, asymptotic dynamics of the perturbation angle 
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ipn in our cases represents a regular two-frequency quasiperiodic motion. 

As the parameters of the maps are varied to a = 0.85, £ = 0.077, the two- frequency torus 
of the map (11) disappears via the backward saddle-node bifurcation, and a three-frequency 
torus appears. A trajectory generated by the map (11) then uniformly fills a 2D torus 
{if E [0,27r),6' G [0,27r)} (see Fig. 4(c)). Strictly speaking, in this situation the map (11) 
does not possess an attractor, since it is characterized by two zero Lyapunov exponents. 
Thus, the dynamics of the perturbation angle is regular and represents a three-frequency 
quasiperiodic motion. 

On the other hand, as the parameters change to a = 0.9, e = 0.28, the two-frequency 
torus of the map (11) destroys via phase-dependent mechanism with the birth of a strange 
nonchaotic attractor shown in Fig. 4(d). Note that the nontrivial Lyapunov exponent re- 
mains negative (cr = —0.7747). The angle variable ipn of an arbitrarily chosen perturbation 
in a vicinity of the torus T asymptotically behaves as (</?„ — $(^n)) ~ exp{an) as n — > oo, 
where $(^) is a fractal- like function characterized by non-differentiability and upper /lower 
semi continuity [4]. As a result, the dynamics of the perturbation angle ipn is irregular. A 
graphic image of such dynamics is a trajectory on a strange nonchaotic attractor of the 
map (11). 

In order to observe a configuration of regions corresponding to different types of dynamics 
of the perturbation angle variable (f in the a — e parameter plane we need to combine the 
results obtained for the attractors of the maps (1) and (11). Such a combined phase diagram 
is presented in Fig. 5(a). The regions of chaos, SNA and chaotic transient in the map (1) are 
shown in black, dark-grey and pattern, respectively (the same as in Fig. 2(a)). The regions 
of the existence of a smooth attractor (i.e. T, 2T and T2 ) of the map (1) are subdivided 
according to the following principle: the white tone (1,4) denotes regions of three-frequency 
quasiperiodic dynamics of the angle variable (p (three-frequency torus in the map (11)), the 
light-grey tone (2,5,7) shows regions of two- frequency quasiperiodicity of the perturbation 
angle ip (two-frequency torus in the map (11)), and the grey tone (3,6,8) corresponds to 
regions of irregular dynamics of the angle <^ (SNA in the map (11)). An enlarged fragment 
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of the diagram in the region of complex dynamical transitions in the model map (1) is shown 
in Fig. 5(b). 

From Fig. 5(a) one can see that the configuration of regions of different behaviour of Lp is 
rather specific. The regions of a two- or three-frequency quasiperiodicity and of a SNA are 
organized as Arnold tongues. Recall that analogous tongues of the localized and extended 
phases were observed for the Harper equation (Fig. 3). However, the shapes of the tongues in 
Fig. 5(a) are distorted compared with those in Fig. 3. This distortion is associated with the 
non-rigorous transition from the Eq. (13) to Eq. (15), where only first terms of the Fourier 
series (14) were taken into account. 

The line D (T — > 2T) on the phase diagrams (Fig. 5(a),(b)) correspond to the usual 
torus-doubling bifurcation in the map (1). The hne L (T — > T2) denotes the torus length- 
doubling bifurcation, which leads to the birth of a double- wrapped torus T2. The lines 
F, I, M and B in Fig. 5(b) denote the same dynamical transitions on the exit from the 
quasiperiodicity region that in Fig. 2(b). 

Proceeding from Figs. 5(a),(b), let us analyze the dynamics of perturbations on the 
tori T and 2T at the threshold of the complex transitions F, I, M and B. On the phase 
diagrams one can see that all the regions of SNA, chaos and chaotic transient adjoin to the 
quasiperiodicity regions shown in the gray tone (and denoted by the numbers 3 and 6), 
corresponding to the existence of irregular dynamics of the perturbation angle variable </?. 
Thus, we can claim that the onset of irregular behavior of the perturbation angle on the 
torus T precedes the following transitions: destruction of a smooth torus T and the birth 
of a SNA via fractalization (route F) and intermittency (route I) scenarios, destruction of a 
smooth torus T and the appearance of a chaotic transient via collision of the smooth torus 
T with a fractal basin boundary (route B). In the same way, the appearance of irregular 
behavior of the perturbation angle <^ on the double torus 2T precedes the destruction of the 
double torus 2T and the birth of a SNA via coUision of 2T with the saddle parent torus 
(route M). It is important to note that all the transitions mentioned have phase-dependent 
mechanisms, as follows from their analysis in terms of rational approximations [21,24,27]. 
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On the other hand, the hnes of regular bifurcations such as torus doubhng (D) and torus 
length-doubhng (L) separate regions shown in hght-gray tone corresponding to the existence 
of two-frequency quasiperiodic dynamics of the perturbation angle variable (p. Note, that 
regular tori bifurcations have phase- independent character (for instance see Ref. [13]). Thus, 
numerical analysis shows that phase-independent bifurcations of tori in the model map (1) 
are accompanied by regular behavior of small perturbations on a torus. 

In our opinion, the connection between the behaviour of perturbations on a torus and 
the types of torus bifurcations is a problem of significant mathematical interest. Recently, 
two of us ( Jalnine and Osbaldestin) suggested a partial solution of this problem for the case 
of torus-doubling bifurcation [35]. We introduced Lyapunov vectors as basis directions of 
contraction for an element of phase volume in a vicinity of the torus and shown that torus- 
doubling can occur only if the dependence of the Lyapunov vectors upon the angle coordinate 
6 on torus is smooth. We have also shown that the appearance of a non-smooth dependence 
of the Lyapunov vectors upon the angle coordinate on torus terminates the line of torus- 
doubling bifurcation on the parameter plane. Note that, after sufficiently large number of 
iterations, an arbitrarily chosen initial perturbation vector of the form 5vq — {5xq, Si/o, 0) 
on the torus T (see Sec. Ill) in a typical case tends to the leading Lyapunov vector, i.e., 
that corresponding to the largest Lyapunov exponent. Obviously, the regular dynamics 
of perturbations on the torus is associated with a smooth dependence of the Lyapunov 
vectors upon the angle coordinate 9 on the torus, while the appearance of a SNA in the 
dynamics of the perturbation angle (/? manifests the onset of a non-differentiable dependence 
of the Lyapunov vector upon the coordinate 9. The latter observation reveals the reason for 
regularity of the dynamics of (p, which accompanies regular torus bifurcations (e.g. doubling 
or length-doubling of a torus). However, the reason for irregularity in the behaviour of (p 
on the threshold of phase-dependent transitions (such as torus fractalization, intermittency 
etc.) is still unrevealed and worthy of further study. 
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VI. CONCLUSION 



In the present work we investigated the dynamics of small perturbations of trajecto- 
ries on an attracting torus in the quasiperiodically forced Henon map. We have shown 
that dynamics of the angle variable of a perturbation may be a two- or three-frequency 
quasiperiodicity, or it may have irregular character. In the last case the graphic image of 
such dynamics is a trajectory on a strange nonchaotic attractor. It was also shown that 
the appearance of irregular character in the approach of trajectories to an attracting torus 
precedes the destruction of such a torus via phase-dependent mechanisms with the birth of 
a strange nonchaotic attractor or of a chaotic transient. 

ACKNOWLEDGMENTS 

The work was supported by the Russian Foundation of Basic Research (Grant No. 03- 
02-16074), CRDF (BRHE REC-006 ANNEX BF4M06, Y2-P-06-16), and by a grant of the 
UK Royal Society. 



15 



REFERENCES 

[1] C. Grebogi, E. Ott, S. Pelikan, J. A. Yorke, Physica D 13, 261 (1984). 

[2] M. Ding, C. Grebogi, E. Ott, Phys. Lett. A 137, 167 (1989). 

[3] B. R. Hunt, E. Ott, Phys. Rev. Lett. 87, 254101 (2001). 

[4] J. Stark, Physica D 109, 163 (1997). 

[5] G. Keller, Fundamenta Mathematicae 151, 139 (1996). 

[6] S.P. Kuznetsov, A.S. Pikovsky, U. Feudel, Phys. Rev. E 51, R1629 (1995). 

[7] A. Pikovsky, U. Feudel, J. Phys. A: Math. Gen. 27, 5209 (1994). 

[8] A. S. Pikovsky, M. A. Zaks, U. Feudel, J. Kurths, Phys. Rev. E 52, 285 (1995). 

[9] U. Feudel, A. Pikovsky, A. Politi, J. Phys. A: Math. Gen. 29, 5297 (1996). 
[10] J. W. Shuai, K. W. Wong, Phys. Rev. E 57, 5332 (1998). 
[11] A. Pikovsky, U. Feudel, CHAOS 5, 253 (1995). 
[12] K. Kaneko, Prog. Theor. Phys. 69, 1806 (1983). 

[13] H. Broer, G. B. Huitema, F. Takens, B. L. J. Braaksma, Mem. Amer. Math. Soc. 83, 1 
(1990). 

[14] P. R. Chastell, P. A. Glendinning, J. Stark, Phys. Lett. A 200, 17 (1995). 

[15] P. Glendinning, Discrete and Continuous Dynamical Systems - Series B Vol. 6, No 4, 
P. 1 (2002). 

[16] H. Daido, Prog. Theor. Phys. 71, 402 (1984). 

[17] U. Feudel, A. S. Pikovsky, J. Kurths, Physica D 88, 176 (1995). 

[18] H. Osinga, J. Wiersig, P. Glendinning, U. Feudel, Int. J. of Bifurcations and Chaos Vol. 
11, No 12, P. 3085 (2001). 

16 



[19] 


S. 


P. Kuznetsov, Phys. Rev. E 65, 066209 (2002). 


[20] 


T. 


Nishikawa, K. Kaneko, Phys. Rev. E 54, 6114 (1996). 


[21] 


S. 


Datta, R. Ramaswamy, A. Prasad, Phys. Rev. E 70, 046203 (2004). 


[22] 


J. 


F. Heagy, S. M. Hammel, Physica D 70, 140 (1994). 


[23] 


A 


Prasad, V. Mehra, R. Ramaswamy Phys. Rev. Lett. 79, 4127 (1997). 


[24] 


S. 


-Y. Kim, W. Lim, E. Ott, Phys. Rev. E 67, 056203 (2003). 


[25] 


T. 


Yalcinkaya, Y.-C. Lai, Phys. Rev. Lett. 77, 5039 (1996). 


[26] 


H 


M. Osinga, U. Feudel, Physica D 141, 54 (2000). 


[27] 


S. 


-Y. Kim, W. Lim, Phys. Lett. A 334, 160 (2005). 


[28] 





Sosnovtseva, U. Feudel, J. Kurths, A. Pikovsky, Phys. Lett. A 218, 255 (1996) 


[29] 


J. 


Heagy, W. L. Ditto, J. Nonlinear Sci. 1, 423 (1991). 


[30] 


J. 


J. Staghano, J.-M. Wersmger, E. E. Slaminka, Physica D 92, 164 (1996). 


[31] 


A 


Witt, U. Feudel, A. S. Pikovsky, Physica D 109, 180 (1997). 


[32] 


A 


Prasad, R. Ramaswamy, I. I. Satija, N. Shah, Phys. Rev. Lett. 83, 4530 (1999) 


L J 


o. 


kj. IN c^i, jr\. L i dodU., rv. rvciiiicibWciiiiy, JT iiy oicd lj X'±(J, ± j . 


[34] 


S. 


p. Kuznetsov, JETP Lett. 39, 113 (1984). 


[35] 


A 


Yu. Jalnine, A. H. Osbaldestin, Phys. Rev. E 71, 016206 (2005). 



17 



FIGURES 





20 




21 




22 



